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Abstract The aim of this work is to investigate rotation profile 
of solar-like stars with magnetic fields. A diffusion coefficient of 
magnetic angular momentum transport is deduced. Rotating stel- 
lar models with different mass are computed under the effect of the 
coefficient. Then rotation profiles are obtained from the theoretical 
stellar models. The total angular momentum of solar model with 
only hydrodynamic instabilities is about 13 times larger than that 
of the Sun at the age of the Sun, and this model can not reproduce 
quasi-solid rotation in the radiative region. However, not only can 
the solar model with magnetic fields reproduce an almost uniform 
rotation in the radiative region, but its total angular momentum is 
consistent with helioseismic result at the level of 3 a at the age 
of the Sun. The rotation of solar-like stars with magnetic fields 
is almost uniform in the radiative region. But there is an obvious 
transition region of angular velocity between the convective core 
and the radiative region of models with 1.2 - 1.5 Mq, where angu- 
lar velocity has a sharp radial change, which is different from the 
rotation profile of the Sun and massive stars with magnetic fields. 
Moreover the changes of the angular velocity in the transition re- 
gion increase with the increasing in the age and mass. 
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1 INTRODUCTION 

Helioseismology has given us detailed internal information about the structure 
and rotation of the Sun: the Sun's rotation is slow in the core and is almost uni- 
form in the radiative region, but the angular velocity has a latitudinal gradient in 
the convective zone (Gough et al. 1996; Schou et al. 119981 Chaplin et al. |19991 >. 
Although the data of inversion of solar-like stars are limited, it has been revealed 



* E-mail: yangwuming@ynao .ac.cn 



2 



Yang, W. and Bi, S. 



that the localized information of stellar interior can be given by asteroseismology 
(G ough & Kosovich ev[T993l Go ugh [19981 Roxburgh et al. [T998l Berthomieu et 
al. 1200 U Basu et al. 120021 Basu l2003l l. Furthermore, it has been shown that lo- 
calized information on the internal rotation profile of a solar-like star can be ob- 
tained from the frequencies of oscillations (Gough & Kosovichev 1 1 9931 Gough 
[19981 Goupil et al. [19961 Lochard et al. 120041 [2005). The information on the 
internal rotation of (3 Cepheid has already been provided by asteroseismology 
(Aerts et al. 2003). And using the data of the Microvariability and Oscillation 
of Star (MOST) satellite, Walker et al. ( f2007T > found that kappal Ceti has a dif- 
ferential rotation profile closely resembling that for the Sun. With ongoing and 
forthcoming space seismic missions: COvection, ROtatin and planetary Transits 
(COROT) (Baglin |2006l l and Kepler (Christensen-Dalsgaard |2U07l l, it is possible 
to extract the information on the internal rotation profile of solar-like stars. 

Moreover, magnetic fields of the active regions on solar surface are believed 
to originate from strong toroidal magnetic fields generated by solar dynamo at 
the base of the convective zone. The understanding of both the stellar magnetic 
activity and the generation of magnetic fields is dependent on the information 
about the interior rotational properties of stars (Thompson et al. 2003 ; Fan 120041 
Charbonneau 2005 ). However the evolution of rotation profile inside stars is 
poorly understood. Therefore it is an important problem to get a global picture 
of the evolution of rotation profile inside stars. 

The influence of rotation on the stellar structure and evolution is stud- 
ied by many investigators (Kippenhaln & Thomas 119701 Endal & Sofia 
[19761 Pinsonneault et al. [198^1 Meynet & Maeder [19971 Huang et al. [20071 . 
Redistribution of angular momentum within the interiors of stars has also been 
considered by many authors (Endal & Sofia |198Tl Chaboyer et al. 119951 Maeder 
& Mevnet l20"00l Palacios et al. 120031 Huang 120041 . These studies show that 
hydrodynamic angular momentum transport processes are unefficient in stars. 
Therefore magnetic angular momentum transport or other mechanisms should be 
considered in rotating stars. In addition, magnetic angular momentum transport 
in massive stars has been investigated by Maeder & Meynet (2003 , 2004, 2005 ). 
The massive stars with magnetic fields rotate almost as a solid body throughout 
the whole star (Maeder & Meynet 2004). Eggenberger et al. (2005 ) and Yang & 
Bi ((2006 ) study the rotation profile of the Sun and show that the quasi-solid rota- 
tion in the Sun can be achieved by considering the effect of the magnetic fields. 
In this paper, we mainly focus on the internal rotation profiles of solar-like stars. 
In Sect. 2 diffusion coefficient of magnetic angular momentum transport is given. 
In Sect. 3 numerical calculation and results are presented. Then, discussion and 
conclusion are made in Sect. 4. 

2 DIFFUSION COEFFICIENT OF MAGNETIC ANGULAR 
MOMENTUM TRANSPORT 

Spruit ( 119991 12002) developed the Tayler-Spruit dynamo, which can generate 
magnetic fields in the radiative region of differentially rotating stars. These fields 
are predominantly azimuthal components, B ~ B^. If magnetic fields exist in 
stars, magnetic angular momentum transport can be described by magnetic in- 
duction and momentum equations. For a constant magnetic diffusivity and shel- 
lular rotation (Zahn 1992), under axisymmetry and only considering Lorentz 
force, the azimuthal components of the induction and momentum equations are 
(Barnes et al. [19991 Yang & Bi l2006b 



dB^ 



+ V( 2 . 2a - V 2 )B = r sin0B p • VO 



(1) 



dt 



r 2 sin 9 
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pr 2 sin 2 6^ = -^B p • V(r sin 9B^) . (2) 
at An 

If the effect of the magnetic diffusivity is to limit the growth of the toroidal field 
after some time, the growth of the instability is halted by dissipative processes 
that operate on a timescale r. Accordingly the second term on the left-hand side 
of Eq. ([TJ may be replaced simply by /t (Barnes et al. 1999). Substituting for 
the second term of Eq. (Q~|i and differentiating the Eq. CO with respect to time, 
one can obtain (Barnes et al. 1999 ) 

For much longer times than the timescale of the instability, one would expect the 
term involving the first time derivative to dominate (Barnes et al. |19991 l, so that 

i^«,sin.B p .vf . (4) 
t at at 

For shellular rotation 0(r, 6) ~ O(r) (Zahn [19921 , then 

B^ ~ rr sin 6B r (5) 
or 

Using Eq. ©, equation (|2]i can be rewritten as 

Pr 29 l ^^Bp-V^sin 2 ^^) 

-S|(- 2 ^f) (6) 

r z or v 47rp r or ' 

The diffusion coefficient for angular momentum transport can thus be obtained 

as 

rBl 

D m = — ^. (7) 

In paper I (Yang & Bi 2006 ) we also got a similar diffusion coefficient, but it 
was only an assumption that the coefficient can be used in the equation of the 
transport of angular momentum. From Eq. (O it can be found that the magnetic 
angular momentum transport approximately obeys the diffusion coefficient D m . 

For a steady equilibrium, the dissipating timescale r has to match the growth 
timescale cr _1 of the instability. Using the growth time scale of magnetic insta- 
bility given by Pitts & Tayler ( 11985b and Spruit d 19991 ), 

<7 =—,uj a = - — __ (8) 



J 2 '~" (47T P )l/2 r ' 



one can get the diffusion coefficient 



D - ^M- 



(9) 



Equation (O can also be rewritten as 



0J A 

where 



D m = r 2 fi( — ) 2 , (10) 
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Equation (TToT > hints that magnetic angular momentum transport is related to 
Alfven waves. 

The distribution of magnetic fields inside a star is poorly known. One of the 
distributions of magnetic fields was given by Spruit d2002b 

< 12 > 

where q — — ^7-, for the case that the effect of thermal diffusion can be 
neglected, namely > Nt, and 

fr 2,/4 <^ 1/ W ,/4 <»> 

for the case 1 with the effect of thermal diffusion. Using the expressions (TTZl i and 
(TT~3T >. one can rewrite Eq. as 

An = rW(-^) 4 (14) 



for the case and 



D ml = 2 1/ Vr!(A)i/2 ( _^_ ) i/2 (15) 

for the case 1. Equations (fl4l and (15[ are consistent with the effective mag- 
netic viscosity defined by Spruit (2002) and Maeder & Meynet d2004| l for the 
radial transport of angular momentum. The expression of ( TBI and (fT~5T > is only 
one of the cases of D m . Braithwaite (2006) validates the Tayler-Spruit dynamo 
scenario, but which is contrary to the findings of Zahn et al. ([2007 ). The rota- 
tion profile of massive stars with magnetic fields was investigated by Maeder & 
Meynet 2004 And the rotation profile of the Sun with magnetic fields was stud- 
ied by Eggenberger et al. 120051 In this work we focus on solar-like stars with 
mass 1.0 - 1.5 M©. 

3 NUMERICAL CALCULATION AND RESULTS 
3.1 Angular momentum transport and loss 

The Yale Rotation Evolution Code (YREC7) is used to construct stellar models 
in its rotating configuration (Pinsonneault et al. 1989; Guenther et al. 11 9921 ) . All 
models are evolved from fully convective pre-main sequence (PMS) to some- 
where near the end of the Main Sequence (MS). The newest OPAL EOS-20053 
(Rogers & Nayfonov [20021 1, OPAL opacity (Iglesias & Rogers [19961 1, and the 
Alexander & Ferguson ( 1994) opacity for low temperature are used. The mod- 
els take into account diffusion of helium and metals, using the prescription of 
Thoul et al. (119941 1. The initial chemical composition of the models is fixed at 
Z a = 0.02, X = 0.706. 

Hydrodynamic instabilities considered in the YREC7 have been presented 
by Pinsonneault et al. ( 119891 ). It is assumed that convection enforces solid-body 
rotation in the convective regions of a star. Therefore the rotational instabilities 
are effective only in the radiative regions. The transport of angular momentum is 
treated as (Endal & Sofia fl978l Pinsonneault et al. [19891 

pr m= f ^ rD ^ (16) 
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Table 1 Model parameters 



model 


mass 


ra 
J 


fk 


JQl 


fu2 




M Q 


10 50 g crrrV 1 








Ml.Oa 


1.0 


1.591 


3.0 


1.0 





Ml. Ob 


1.0 


1.591 


3.0 


1.0 


0.01 


M1.2 


1.2 


1.9095 


3.0 


1.0 


0.01 


M1.4 


1.4 


1.534 


1.0 


1.0 


0.01 


M1.5 


1.5 


1.095 


1.0 


1.0 


0.01 



a The initial angular momentum; 

b The v alue of fn for the coefficient of hydrodynamic instabilities (Pinsonneault et al. 
[79891 ; 

c The value of fa for D m ; model Ml.Oa with only hydrodynamic instabilities. 



in the radiative regions of a star, where fn is an adjustable parameter introduced 
to represent some inherent uncertainties in the diffusion equation. In stars with 
M < 1.5M Q , angular momentum loss due to magnetic braking is treated as a 
parameterized formula (Kawaler l 1988b 

to reproduce the Skumanich relationship (Skumanich 1972), where Kq ~ 1.13 x 
10 47 g cm 2 s, fx is an adjustable parameter. It is assumed that the magnetic 
braking has an effect on the whole convective envelope. In some case, a PMS star 
could be locked by the surrounding accretion disk. However the disk can extract 
angular momentum from the star as well as can supply angular momentum to 
the star (Stassun & Terndrup 2003). Moreover, it has been argued by Matt & 
Pudritz (2005a 2005b) that the spin-down of PMS stars may not be due to a 
magnetic star-disk interaction, but may result from a magnetic stellar wind. Thus 
for simplicity, we do not consider the magnetic star-disk interaction. 

The initial angular momentum of a star is still uncertain. Kawaler d 19871 ) 
shows that the angular momentum J of stars more massive than 1.5 Mq is pro- 
portional to squared mass M 2 . But the mass-momentum relation of stars below 
1 .5 M Q is uncertain. As a first test, we take the initial angular momentum to be 
a free parameter. The adjustable parameters mentioned above are listed in Table 

m 

3.2 Results of calculation 

FigureQlcompares the evolution of the internal rotation profile of 1.0 M mod- 
els with only hydrodynamic instabilities given by Pinsonneault et al. ( 1989]) and 
with both the hydrodynamic instabilities and magnetic fields. Both models are 
evolved from PMS with initial angular momentum Jo = 1.591 x 10 50 g cm 2 
s _1 to the age of 4.5 Gyr. During the PMS phase, although with angular momen- 
tum loss from the surface of models, the rotation rate rapidly increases due to 
quickly contracting. The internal rotation profile of model with magnetic fields 
has been different from that of model with only hydrodynamic instabilities when 
the models are near the zero-age main sequence (ZAMS). The model with only 
hydrodynamic instabilities has a fast rotation core. The rotation of model with 
magnetic fields is, however, almost uniform. During the early stage of the MS, 
the rotation of the model Ml.Oa is differential, but the model Ml. 0b is a quasi- 
solid body rotation. At the age of 4.5 Gyr, the surface rotation rates of two models 
are around 2.7 x 10~ 6 rad/s. However, internal distribution of angular velocity is 
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Fig. 1 Rotation profiles as a function of radius for 1 .0 M Q model at 
different ages labeled by Gyr. The initial angular momentum Jo = 
1.591 x 10 50 g cm 2 s -1 . The dotted line indicates the solid-body ro- 
tation of the models with and without magnetic fields on the pre-main 
sequence. 



quite different. The model Ml.Oa shows a strong differential rotation with a fac- 
tor of about 40 between the angular velocity in the core and at the surface; but the 
model Ml. 0b shows an almost uniform angular velocity, with a small increase in 
f2(r) in the center of r < 0.2 Rq, as that obtained by Eggenberger et al. (2005 ). 
The surface rotation rate of model Ml. 0b is higher than that of model Ml.Oa in 
the early evolutionary stage. And the loss rate of angular momentum is related to 
f2 3 . Consequently the amount of the angular momentum loss of model Ml. 0b is 
larger than that of model Ml.Oa. The total angular momentum of model Ml.Oa 
is 2.628 x 10 49 g cm 2 s _1 at the age of 4.5 Gyr, which is about 13 times larger 
than the seismical result (1.94 ± 0.05) x 10 48 g cm 2 s _1 (Komm et al. 120031 ); 
but the total angular momentum of model Ml. 0b is 2.045 x 10 48 g cm 2 s — 1 at 
the age of 4.5 Gyr, which is consistent with the result of helioseismology at the 
level of 3 a. 

Figure [2] shows the evolution of the internal rotation profile of model Ml. 2. 
In the early evolutionary stage, the angular velocity f2(r) is almost constant. At 
the stage of X c ~ 0.69, about 90 percent of the initial angular momentum has 
been lost; the rotation is nearly uniform in the radiative region; but the rotation of 
the convective core is faster than that of the radiative region. There is a transition 
region between the convective core and the radiative region, where the angular 
velocity has a sharp radial change due to the spin-down of the outer parts of 
model resulting from angular momentum loss and the decrease of the horizontal 
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— 1 .227 
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0.706 — 1.829 


701- 


-0.407 


\ 0.690- 


-0.176 
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0.058- 


-0.039 



0.5 1 1.5 2 2.5 



Fig. 2 Rotation profiles as a function of radius for model Ml. 2 at the 
different evolutionary stages indicated by the central hydrogen X c . The 
initial angular momentum Jo = 1.90 95 x 10 50 g cm 2 s _1 . The J to t is the 
total angular momentum of models in 10 50 g cm 2 s _1 . The lower dotted 
line shows the rotation profile of PMS model at the age of 1 Myr. 
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Fig. 3 Same as Fig.[2]but for model Ml .4. The initial angular momentum 

Jo = 1.5 34 x 10 50 gem 2 s" 1 . 



coupling provided by the magnetic field resulting from the increase in /i-gradient 
in the region. The loss of angular momentum mainly occurs in the early evolu- 
tionary stage. During the late stage of the MS, the total angular momentum of 
model Ml. 2 is only several percent of the initial angular momentum; the rota- 
tion is slow; thus the loss rate of angular momentum is very low. Consequently, 
the angular momentum of model Ml. 2 is almost conservative from the stage of 
X c = 0.153 to the stage of X c = 0.058. 

The evolution of the internal rotation profiles of models Ml. 4 and Ml. 5 are 
shown in Figs. [3] and [4] The model Ml. 4 has lost about 50 percent of the ini- 
tial angular momentum at the stage of X c — 0.675. But the model Ml. 5 has 
only lost about 15 percent of the initial angular momentum even at the stage of 
X c = 0.487. The distributions of the angular velocity of models Ml. 4 and Ml. 5 
are different from that of model Ml. 2. The angular velocity between the convec- 
tive core and the radiative region of models Ml. 4 and Ml. 5 decreases obviously 
when the radius increases. The radial change of the angular velocity between the 
convective core and the radiative region in models Ml. 4 and Ml. 5 is larger than 
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Fig. 4 Same as Fig.|2]but for model Ml .5. The initial angular momentum 

Jo = 1.095 x 10 50 gem 2 s" 1 . 




Fig. 5 Internal distribution of the hydrogen mass fraction X as a function 
of the radius. The solid lines show model Ml. 2. The dotted lines refer to 
model Ml. 4. 



that in model Ml .2, which should be due to the /i-gradient and the fast spin-down 
occurring at the same stage in models Ml. 4 and Ml. 5. However, in model Ml. 2, 
the fast spin-down occurs in the early evolutionary stage when the /i-gradient is 
small. 

Figure [5] shows the distribution of the hydrogen mass fraction X of models 
Ml. 2 and Ml. 4. It is obvious that there is a sharp /i-gradient at the bottom of 
the radiative region of models Ml. 2 and Ml. 4 at the late stage of MS. The /i- 
gradient and the O-gradient are in the same region. The ratio of magnetic field, 
B r /B,p, is related to V^ 1 in Tayler-Sprut dynamo, namely D rn ~ V~ 2 . Thus 
the increase in /i-gradient must lead to the decrease in the coupling provided by 
magnetic fields. This scenario was first found by Eggenberger et al. (2005 ) in the 
Sun. 

4 DISCUSSION AND CONCLUSIONS 

The surface velocity is sensitive to the loss rate of the angular momentum. The 
equation ( TP71 ) may overestimate the loss rate of angular momentum of the rapid 
rotation stars (Andronov et al. 120031 1. But it can reproduce the Sun's rotation. 
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Thus we take it in our models. The adjustable parameter, fk, is adjusted to obtain 
the solar rotation rate at the age of 4.5 Gyr in models Ml .Oa and Ml .Ob. However 
the same value of cannot be applied to models Ml. 4 and Ml. 5 because the 
convective envelope of models Ml. 4 and Ml. 5 is too shallow. Thus we take a 
small f k for models Ml. 4 and Ml. 5. 

The value of parameter /n2 is adjusted to obtain a quasi-solid rotation in our 
models. The value of 0.01 can do work in our models. However this value is less 
than 1 . This could be a consequence of overestimating the ratio of B r to B^, in 
Tayler-Spruit dynamo. 

The distribution of angular velocity of models Ml.Oa and Ml. Ob shows that 
the rotation profile strongly depends on the efficiency of angular momentum 
transport. The angular momentum is effectively transported outward by mag- 
netic fields in Ml. Ob. Thus the rotation of the core of Ml. Ob is slow comparing 
with that of Ml.Oa. The surface rotation rate mainly depends on the loss rate of 
angular momentum and the amount of outward transport of angular momentum. 
Because both Ml.Oa and Ml. Ob have same value of fk and initial f2, the discrep- 
ancy of the surface velocity between Ml .0a and Ml .Ob relies on the efficiency of 
outward transport of angular momentum. In Ml. Ob, the loss of angular momen- 
tum is counteracted by magnetic angular momentum transport. Thus the surface 
velocity of Ml. Ob is higher than that of Ml.Oa when the interior of Ml. Ob has 
enough angular momentum to transport outward. The loss rate of angular mo- 
mentum is related to fi 3 . Consequently the amount of angular momentum loss 
of model Ml. Ob is larger than that of model Ml.Oa. This scenario takes place in 
the early evolutionary stage. 

At the early stage of Ml .2, the fast spin-down leads to the sharp radial change 
of angular velocity at the top of the convective core. But at the same stage of 
models Ml. 4 and Ml. 5, the spin-down is slow. At the late evolutionary stage of 
Ml. 2, although there is a large /i-gradient at the top of the core, the spin-down 
is very slow. Thus the radial change of the angular velocity is small comparing 
with that of models Ml. 4 and Ml. 5 at the top of the core. However in models 
Ml. 4 and Ml. 5, the /i-gradient and the fast spin-down resulting from angular 
momentum loss and stellar expansion occur on the same stage. Therefore the 
radial change of angular velocity is large at the top of the core in models Ml. 4 
and Ml. 5. 

The 1 .0 M Q model with only hydrodynamic instabilities has a fast rotation 
core, and its total angular momentum is 2.628 x 10 49 g cm 2 s -1 at the age of 4.5 
Gyr, which disagrees with the helioseismic results. However the 1.0 M Q model 
with magnetic fields has a slow rotation core, and the rotation is almost uniform 
in the radiative region, which are consistent with the seismical results. Moreover 
the total angular momentum of the model with magnetic fields is 2.045 x 10 
g cm 2 s~ 1 at the age of 4.5 Gyr, which agrees with the helioseismic result at the 
level of 3 a. 

A diffusion coefficient of magnetic angular momentum transport is obtained. 
Not only can the magnetic fields reproduce a quasi-solid rotation, but they can 
enhance the loss rate of angular momentum. The rotation of solar-like stars with 
magnetic fields is almost uniform in the radiative regions, which is consistent 
with the results of helio- and asteroseismology. However there is a transition 
region between the convective core and the radiative region, where the angular 
velocity has a sharp radial change, which is different from that of solar model and 
that of massive stars shown by Maeder & Meynet (2004 ). Moreover the changes 
of the angular velocity in the transition region increase with the increasing in the 
age and mass. 
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